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Experiments by Gittings, Bandyopadhyay, and Durian [Europhys. Lett. 65, 414 (2004)] demon- 
strate that light possesses a higher probability to propagate in the liquid phase of a foam due 
to total reflection. The authors term this observation photon channelling which we investigate in 
this article theoretically. We first derive a central relation in the work of Gitting et al. without 
any free parameters. It links the photon's path-length fraction / in the liquid phase to the liquid 
fraction e. We then construct two-dimensional Voronoi foams, replace the cell edges by channels 
to represent the liquid films and simulate photon paths according to the laws of ray optics using 
transmission and reflection coefficients from Fresnel's formulas. In an exact honeycomb foam, the 
photons show superdiffusive behavior. It becomes diffusive as soon as disorder is introduced into the 
foams. The dependence of the diffusion constant on channel width and refractive index is explained 
by a one-dimensional random-walk model. It contains a photon channelling state that is crucial for 
the understanding of the numerical results. At the end, we shortly comment on the observation that 
photon channelling only occurs in a finite range of e. 

PACS numbers: 05.40.Fb,82.70.Rr,42.68.Ay,05.60.-k 



I. INTRODUCTION 

Aqueous foams consist of gas bubbles separated by 
liquid films 0. In dry foams, these bubbles are de- 
formed to polyhedra. Always three films meet in the 
so-called Plateau borders which form a network of liq- 
uid channels throughout the foams. Their opaque ap- 
pearance identifies them as multiply scattering media. 
Moreover, careful light-scattering experiments show that 
light transport has reached its diffusive limit in foams 
Qi ^ lol which means that photons perform a ran- 
dom walk. Diffusion of light is well established in col- 
loidal systems , where light is scattered from particles, 
or in nematic liquid crystals 8], where scattering occurs 
from fluctuations in the local optical axis. Diffusing- 
wave y and diffuse-transmission spectroscopy pro- 
vide non-invasive probes of structure and dynamics in 
bulk samples. While these methods have already suc- 
cessfully been applied to monitor the internal dynamics 
of foams Q, H, Q 1^ , no clear understanding has emerged 
so far about the basic mechanism underlying the random 
walk of photons in these cellular structures. Estimates 
for the scattering from Plateau borders do not seem to be 
in accordance with experiments ^| . Scattering from 
vertices formed by four meeting Plateau borders would 
be another possibility On the other hand, since gas 
bubbles in foams are much larger than the wavelength 
of light, one can use the laws of ray optics to follow the 
photons on their random walk as they are reflected by 
the liquid films jll [H El El . 

The work presented in this paper was initiated by a 
very inspiring publication of Gittings, Bandyopadhyay, 
and Durian published under the same title The 
authors argued that total reflection of light at the liquid- 
gas interface in a foam should increase the probability of 
photons to propagate in the liquid phase, made up by 
thin films and Plateau borders, and they termed this ef- 



fect photon channelling. By measuring the absorption of 
light in foams with strongly absorbing liquid, they deter- 
mined the fraction / of a photon's path that lies in the 
liquid. Equal distribution of the photons throughout the 
foam would give / = e, where e is the liquid fraction in 
the foam. However, the authors found that in a certain 
range of e the path-length fraction / fulfills f > e. This 
not only confirms photon channelling but it also gives 
clear evidence that the basic mechanism for light trans- 
port in foams depends on the liquid volume fraction. 

Motivated by these findings, we perform a thorough 
theoretical study of photon channelling in this article 
based on analytical arguments and numerical simula- 
tions using ray optics. Gittings et al. derived a rela- 
tion between / and e E3- It contains the ratio of av- 
eraged transmission probabilities which they determined 
both by experiments and simulations. In Sec. IIII Al we 
demonstrate how this ratio follows from general argu- 
ments and are therefore able to justify the relation /(e) 
on completely analytical grounds. Furthermore, Gittings 
et al. stated that the effect of photon channelling would 
prevent the photons from performing a "truely random 
walk" Uoj. Indeed, our numerical investigation in Sec. 
IIV Al shows that transport of light in a perfect honey- 
comb foam is superdiffusive. The effect is most pro- 
nounced when the liquid-gas interface is totally reflecting 
|l6j |. However, the transport becomes diffusive as soon as 
the foam contains disorder. This demonstrates that pho- 
ton channelling is compatible with the concept of ran- 
dom walk. To complete our numerical study of photon 
channelling, we have developed a one-dimensional ran- 
dom walk model. It contains a photon-channelling state 
and can therefore explain the main features of our nu- 
merical investigation. 

The paper is organized as follows. In Sec.m we present 
two-dimensional model foams for our photon-channelling 
study and explain how we simulate the photons' random 
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walk based on the laws of ray optics. In Sec. lIIII charac- 
teristic features of photon channelling are first explained 
analytically and then compared with simulation results. 
The transport properties of light in our model foams are 
described in Sec. II VI Numerical results on supcrdiffusion 
in perfect honeycomb structures and diffusion in disor- 
dered foams are presented. Finally, we introduce the one- 
dimensional random- walk model for photon channelling. 
Section Ivl summarizes and discusses our results. 



II. MODEL SYSTEM AND METHOD OF 
SIMULATION 

To introduce a simple model for a two-dimensional 
foam, we start from a Voronoi tessellation of the plane 
[l7lll8l |. It is generated from a distribution of seed points 
in a simulation box, for which Voronoi polygons are con- 
structed in complete analogy to the Wigner-Seitz cells 
of periodically arranged lattice sites. We begin with 
a triangular lattice of seed points that gives a regular 
honeycomb structure whose edges possess the length /q. 
Then we systematically introduce disorder by shifting the 
seed points along a randomly chosen displacement vec- 
tor whose magnitude is equally distributed in the inter- 
vall [0,(5r]. Examples of Voronoi foams are presented, 
e.g., in Ref. jl3|. Note, while the structure and evolu- 
tion of real foams are constrained by Plateau's rules and 
Laplace's law 0, simple honeycomb, Voronoi, and three- 
dimensional Kelvin structures were used for a first study 
of the physical properties of foams 19, 20]. 

All our Voronoi tessellations are produced by the soft- 
ware Triangle They contain approximately 15000 
cells which corresponds to a quadratic simulation box 
with edge length 200 Iq. Now, on each edge of the Voronoi 
tessellation, we place a channel of width d that represents 
a liquid film. Our final model foam for (5r = is indicated 
in Fig.^ For this regular honeycomb foam, we calculate 
the fraction of the plane filled by the liquid phase as 

In the following, we only investigate foams with modest 
disorder quantified by Sr < O.SZq. This avoids the un- 
physical situation that four instead of three channels or 
films meet when we construct our model foams. All cells 
still have six edges. Furthermore, to simulate photon dif- 
fusion in the model foams, periodic boundary conditions 
are used. 

We let photons perform a random walk in our model 
foams by applying the rules of geometrical optics (see Fig. 

The photons move straight in the liquid or gaseous 
phase with respective velocities c/ni or c/ng where c is 
the vacuum speed of light and ni , ng are refractive in- 
dices. When the photons hit a liquid-gas interface, they 
are reflected with a probability, also called reflectance, 
given by Fresnel's formulas. The respective reflectances 




FIG. 1: Regular honeycomb foam {5r = 0) and a photon's 
random walk in this structure. 

for electric polarizations parallel or perpendicular to the 
plane of incidence are 

^ / tan(A-/3,) y ^ / sin(A-/?,) y 

II \t&n{(3i + f3g) J ^ \smiPi + (3g) J ' 

(2) 

where /?; and (3g denote the respective angles of a light 
ray in the liquid or gaseous phase measured relative to 
the normal on the interface. They are connected by Snel- 
lius' law: sin/?/ /sin/3g = Ug/rii. Note that Fresnel's for- 
mulas are symmetric in /?; and Pg. So the reflectance is 
the same whether the light ray is hitting the interface 
coming from the liquid or the gaseous phase. However, 
there is an important difference: light rays in the opti- 
cally denser medium (n; > ng) experience total reflection 
for an incident angle Pi > f3* = arcsm{n g/ni). We stress 
this point here because it is central for the occurence of 
photon channelling. 

Typically, we launch 10000 photons at one vertex of 
the underlying Voronoi tcsselation in an angular range 
of 60° and let them run during a time t — IQ^Iq/c. At 
several times, we calculate the mean square displacement 
(7^ from the photon cloud and plot it as a function of 
t. Diffusion constants are then determined from a flt to 
= iDt. 



III. PHOTON CHANNELLING 

At the core of the observation of photon channelling 
is the fact that photons can totally be reflected when 
they hit the liquid-gas interface from the liquid side. Ac- 



3 



cordingly, a photon spends more time in the Hquid phase 
or, more accurately, the fraction / of the photon's path 
that hes in the Uquid phase is not equal to the Uquid 
volume fraction e. Gittings et al measured / with the 
help of a strongly absorbing liquid and then made a 
prediction for the two averaged transmission probabil- 
ities for photons going from the liquid to the gaseous 
phase (Ti^g) or the reverse direction (Tg^i). Only in 
the range 0.04 < e < 0.2, they found that the ratio 
Ti^g/Tg^i deviates from one, indicating photon chan- 
nelling. Furthermore they gave its approximate value 
as Ti—,g/Tg—,i ^ Ug/ni and also confirmed it by simula- 
tions. In the following, we derive this central relation for 
Ti^g/Tg^i analytically, connect it to / and compare it 
with our simulations. 



A. Analytic results 

We start with defining the averaged transmission prob- 
abihty Ti^g by 



tt/2 



Pl{Pl)t{Pl,Pg)d/3l 



(3) 



where the angle-dependent transmission coefficient 
t{Pi,f3g) = I — r{Pi, f3g) is connected to the reflectance 
r of Fresnel's formulas introduced in Eq. (|2Jl. Note that 
only the angular range of f3i < f3* contributes to the inte- 
gral since t{f3i,Pg) is for (3i > /?* due to total reflection. 
The factor pi{pi) is the normalised probability of a pho- 
ton to reach the interface from the liquid side with an 
incident angle f3i. We write it as 



PiiPi) = PiiPi)/ 



7t/2 



Pi{fii)dl3i 



(4) 



Accordingly, we define the transmission probability Tg^i 
by 



Pg{Pl)t{(3uPg)dPl 



(5) 



where we used the symmetry of Fresnel's formulas with 
respect to /3; ^ (3g. We still take the angle fii for the 
averaging, therefore the upper limit of the integral is just 
(3* since photons entering the liquid cannot have a /3; 
larger than (3* . For (3i < (3* it makes sense to assume 
that photons coming from the gaseous into the liquid 
phase possess the same angular distribution as photons 
in the liquid, i.e., Pg{Pi) oc Pi{Pi). So we write 

PgiPl) = Pg{Pl)l Pg{(3l)d(3l 

Ja 

= Pimi / Piimpi ■ (6) 

Jo 

Note again the different upper limit in the integral of the 
normalization factor compared to Eq. 10}. Since t{f]i, (3g) 



is for (3i > P* , the transmission probabilities of Eqs. 
Q and (0 appear the same however they differ by the 
normalization factors of pi and pg in Eqs. Q and © 
as just mentioned. Therefore the ratio of the averaged 
transmission probabilities is 



Pi{f3,)d(3i 



T, 



(7) 



We now assume that all directions of the photons are 
equally distributed in the liquid and therefore obtain 
for the angular distribution of photons reaching the in- 
terface, Pi{Pi) = cos Pi. The term cos/?; means that 
the number of photons hitting a given surface area de- 
creases with increasing i.e., photons approaching the 
interface on a "shallow" path {(3i 7''/2) have a van- 
ishing probability to hit a given surface area. With 
/3* — arcsin(ng/n;), the ratio of the mean transmission 
probabilities becomes immediately 



Ta^l 



COS Pid/3i _ n 
j^"^^^ cos /3id(3i 



g 
ni 



(8) 



This is the result of Gittings et al. [T3| measured in the 
range 0.04 < e < 0.2 for the liquid fraction. In our simu- 
lations we realize Pi {f3i ) = cos Pi by starting the photons 
in the liquid with a uniform angular distribution. In ex- 
periments, photons always enter from the gaseous phase. 
A fast randomization of the angular distribution in the 
liquid phase is then achieved by disorder in the foam and 
by the strong curvature of the interface at the Plateau 
borders. Such a randomization due to curvature occurs, 
e.g., in the Sinai billard [2.2j. 

We now derive how the fraction / of a photon's path 
that lies in the liquid phase depends on the liquid volume 
fraction e. We assume a stationnary photon distribution 
in the foam. Then the photon current densities ji^g and 
jg—,1 across an interface have to be equal: ji^g = jg—>i- 
With ji^g = {c/ni)QiTi.^g and jg^i = {c/ng)ggTg^i, 
where Qi is the number density of the photons in phase 
i, and using relation (jSJ, we find immediately 



Qg 



(9) 



The ratio of the average times ti and Tg a photon spends 
in the liquid or gas equals the ratio of the average number 
of photons in these phases. 



QlVl 
Qg^g 



(10) 



where Vi and Vg denote the respective volumes. To moti- 
vate this equation, we consider the average photon num- 
ber as a time average over the paths of many photons 
so that Ti cx QiYi. In experiments, the average path- 
length fraction / — si/{si + Sg) is measured. Here Si is 
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the path length in phase i. With Si = {c/ni)ti and Eqs. 
©, lUni), we finally obtain 



Sl 



n„ / 1 
1 + ^ 



Sl 



ni \e 



(11) 



where e = Vz/(Vz + V^). With the actual values of Ug = 1 
and ni — 1.33, i.e., Ug/ni = 0.75, we obtain the semi- 
empirical relation for /^^ that Gittings et al. report in 
Ref. TcH. 

If we do not specify the value of Ti^g/Tg^i in the 
derivation of Eq. (|ll|l . we obtain the relation 



Ti^g ^ g(l - /) 
Tg^i /(I - e) 



(12) 



that was already given by Gittings et al. It links 

Ti^g/Tg^i to the measurable quantity /. 
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FIG. 2: a) Inverse of the path- length fraction / as a function 
of liquid fraction e for a disordered model foam {Sr/lo = 0.3) 
and a perfect honeycomb structure {Sr/lo ~ 0.0). Photons 
with parallel or perpendicular polarization states are chosen. 
The full line corresponds to Eq. 111^ and the dashed line to 
— e~^. Further parameters are Ug — 1 and ni — 1.33. b) 
The ratio of the average transmission probabilities Ti^g/Tg^i 
as a function of e. Otherwise, same description as in a). 



B. Numeric results 

We performed simulations to confirm the analytic re- 
sults of the previous section always choosing Ug = 1. 
For ni = 1.33, we plot in Fig. EJi) the inverse of the 
path-length fraction / as a function of the liquid frac- 
tion e which was calculated via Eq. from the channel 
width d/lo. A perfectly ordered model foam and one with 
disorder 5r = 0.3lo were chosen. The photons were ei- 
ther in the parallel or perpendicular polarization state. 
The numerical results represented by the symbols agree 
very well with the analytic prediction of Eq. (|lll) (full 
line). Therefore, the diagram and also further simula- 
tions clearly show that the validity of Eq. (|ll|l is indepen- 
dent of disorder in the model foam and the polarization 
state. The dashed line corresponds to / = e. Accord- 
ing to Eq. 111|) this means ng = ni, i.e., when there is 
no optical contrast between the cells and the channels. 
Differently speaking, photon channelling should always 
occur in this model for Ug ^ ni. In Fig. 03), we plot 
Ti^g/Tg-,1, calculated from Eq. (|12|l . as a function of 
£ using the results of Fig. |2t). Again the different nu- 
merical data points agree with our theory culminating in 
Eq. ((SJ and therefore justify the assumptions made dur- 
ing the derivation of Eq. |(SJ). Finally, we also checked 
that the effect of photon channelling becomes stronger 
for increasing n; as expected. 

IV. PHOTON DIFFUSION 
A. Superdiffusion in Honeycomb Foams 

Here we investigate the behavior of photon propaga- 
tion in the exact honeycomb foam of Fig. ^ In Fig. 
we plot the mean-square displacement of the photon 
cloud (in units of Iq) as a function of time t (in units 
of Iq/c) for different refractive indices ni of the liquid 



phase. The channel- width-to-length ratio is d/lo — 1. 
As a reference, the dashed line indicates pure diffusive 
behavior with oc t. Clearly, in the exact honeycomb 
foam the spreading of the photons is superdiffusive. The 
effect is strong for short times with an exponent m in 
a'^ (X between 1.5 and 1.6. It becomes weaker for 
large times, e.g., m = 1.1 for n; = 1.33. In the frame- 
work of Levy walks, superdiffusion is associated with a 
distribution of step lengths whose first or second moment 
(mean value and variance) do not exist When the 

boundaries of the liquid channels are completely reflect- 
ing, photon paths occur that correspond to effective steps 
of infinite length along the directions of the channels [Tg ■ 
In such a system, superdiffusion is most pronounced also 
for large times as we shall demonstrate in a forthcom- 
ing publication 0|. It becomes weaker by letting the 
photons enter the gaseous phase since this is a mean to 
reduce the number of very long or infinite photon steps. 
Note that in Fig. 13 superdiffusion is less pronounced for 
smaller n/, i.e., when the reflectance at the liquid-gas in- 
terface decreases. Another mean to cut long steps would 
be rounding off the sharp edges in the channel system. 
Finally, as demonstrated in the next section, introducing 
disorder in the model foam leads to conventional photon 
diffusion as observed in experiments. 

In Fig.0I we plot the temporal evolution of the mean- 
square displacement for different ratios d/lQ of channel 
width to length, the refractive index is ni = 2. For large 
ratios such as d/lo = 1, superdiffusion is stronger since 
then photons have more possibilities to realize very long 
effective steps. These steps are visible in the insets that 
show the photon clouds at time t = IO^Iq/c for d/lo = 1 
and 10^"^ (see also Fig. EJ. The clouds are not spheres 
as for conventional diffusion, instead, due to the long 
effective steps, they exhibit peaks along the six equivalent 
directions of the honeycomb foam. The peaks are more 
pronounced for stronger superdiffusion. 
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FIG. 3: Exact honeycomb foam: Reduced mean-square dis- 
placement cr'^/io of the photon cloud as a function of reduced 
time t/{c/lo) for different refractive indices n; of the liquid 
phase and for d/lo = 1. Perpendicular polarization is cho- 
sen. The dashed line indicates pure diffusion. Insets: Photon 
clouds at t — lO^lo/c for n; = 1.33 (top left) and 3.0 (bottom 
right). 
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FIG. 5: Diffusion constant D in units of Iqc as a function of 
disorder 5r/lo in the model foam with d/lo ~ 0.1, Ug — 1 and 
71; = 1.33. 
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the logarithmic scale for Sr) that the divergence is quite 
weak especially for the perpendicular polarization state. 
On the whole range of Sr, a clear decrease of the dif- 
fusion constant with increasing disorder is visible. We 
understand this observation since correlations between 
successive photon steps are destroyed more easily in a 
disordered foam. Furthermore, the diffusion constant 
for parallel polarization is larger compared to the per- 
pendicular case due to its smaller reflectance r|| , which 
even becomes zero at the Brewster angle characterized 
by Pi+Pg = n/2. 

In Fig. the dependence of the diffusion constant on 
the channel width to length ratio d/lo is shown for the 
most disordered foam in our numerical treatment, i.e., 
for Sr/lo = 0.3. The ratio d/lo is directly connected to 
the liquid fraction e via Eq. which offers an expla- 
nation for the observed decrease of D with increasing 
d. When the liquid fraction s becomes larger, the pho- 



FIG. 4: Exact honeycomb foam: Reduced mean-square dis- 
placement cr^/Zo of the photon cloud as a function of re- 
duced time t/{c/lo) for different ratios d/lo. The refractive 
index is n = 2 and perpendicular polarization is chosen. The 
dashed line indicates pure diffusion. Insets: Photon clouds at 
t = lO^lo/c for d/lo = 1 (top left) and 10"^ (bottom right). 



B. Diffusion in Disordered Voronoi Foams 

As soon as we introduce disorder into the Voronoi 
foam, superdiffusion vanishes and conventional diffusion 
occurs. This is not unexpected since disorder inhibits 
long effective steps of a photon as explained in the pre- 
vious section. In Fig. [S] we plot the diffusion constant D, 
for both parallel and perpendicular polarization of the 
light wave, as a function of 6r, which quantifies disorder 
in our Voronoi foams. We expect the diffusion constant 
to diverge when Sv/Iq because of the superdiffusiv 
behaviour at 6r/lo — 0. Our simulations indicate (note 
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FIG. 6: Diffusion constant D in units otloc as a function of the 
channel width to length ratio d/lo in a disordered model foam 
with 5r/lo — 0.3, Ug = 1 and n; — 1.33. Symbols indicate 
numerical results, the lines follow from Eq. 1151 derived by a 
a one- dimensional analytical model. 
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FIG. 7; Diffusion constant D in units of Iqc as a function of 
refractive index ni of the liquid phase in a disordered model 
foam with 5r /Iq — 0.3, d/lo = 0.1 and n; = 1. Symbols indi- 
cate numerical results, the lines follow from Eq. 1151 derived 
by a one-dimensional analytical model. 

tons spend more time in the liquid phase where they 
move with a smaller velocity compared to the gaseous 
phase. Furthermore, more photons exhibit total reflec- 
tion at the liquid-gas interface, i.e., photon channelling 
is more pronounced. Note also that for d, e ^ 0, the dif- 
fusion constant approaches a finite value. This is in clear 
contrast to experiments and we will comment on it in 
our conclusions. We have developed a one-dimcnsional 
random-walk model which includes a photon-channelling 
state and whose details are explained in the following sec- 
tion. Its prediction for the diffusion constant (see lines in 
Fig. gives a remarkable quantitative agreement with 
the numeric results in Fig. El Finally, Fig.[7|demonstrates 
that the diffusion constant decreases with increasing n;. 
This agrees with intuition since the reflectance increases 
and the photons in the liquid phase move slower. The 
random-walk model (see lines in Fig. [7|) again gives a 
good description for the observed behavior with some 
quantitative deviations for the parallel polarization state. 



C. A One-Dimensional Model for Photon 
Channelling 

To explain the main features of photon channelling in 
disordered foams, we have developed a one-dimensional 
random-walk model. We consider a sequence of gaseous 
and liquid regions with a lattice constant Ax, as illus- 
trated in Fig. |S1 On a length scale larger than Ax, we 
then introduce the coarse-grained probabilities p'^{x,t), 
p'g'{x, t), Pl*{x, t) and p^{x, t) for a photon to be at time 
t at the position x and in a state specified by the indices. 
Here the lower index means that the photon resides in the 
gaseous (g) or liquid (I) phase and the arrow of the up- 
per index gives the direction of motion. All photons are 
transmitted through the liquid-gas interface with a prob- 
ability T or they are reflected with a probability 1 — T. 
In Fig.|Slall possible transitions between the four states 
are indicated by solid lines, the transmission probabilities 
are given in non-italic writing. Reflected photons return 
to their original position, which for concreteness we lo- 
cate in the middle of each region, after a time Tg or t;. 
Transmitted photons move from the gaseous to the liquid 
phase or vice versa during the time Tgi — {Tg + t;)/2. 

So far, all photons can traverse the liquid-gas interface 
with a probability T. To allow for photon channelling in 
our model, i.e., photons which exhibit total reflection, we 
introduce a fifth state located in the liquid phase with the 
probability pi^dx^t). Within the liquid phase, photons 
travelling to the right or left switch instantaneously into 
the photon channelling state with a probability a. They 
return to one of the "moving" states with a probability 
b/2. The numbers a and b are related to each other as 
we will demonstrate below. With probability 1 — 6, a 
photon in the photon channelling state is reflected at the 
liquid-gas interface and returns to the same state after 
the time t; . All possible transitions connected to photon 
channelling are indicated by dashed lines in Fig. |S1 

The following master equations now quantify all tran- 
sitions between different photon states: 





= {l-T)p^{x,t- 


-rg)^ 


-T{1 - 


- a)p'l* [x — Ax, t — 








= il-T)p^{xA- 




-T{1- 


- a)pp(a; + Ax, t — 




- ^bpiAx,t) 


pT{x,t) 


= Tp'g*{x — Ax, t - 


- Tgi) 




r)(l - a)pf{x,t- 




pr{x,t) 


— Tp'g'{x + Ax,t - 


-Tgi) 


^(1- 


r)(l - a)pi^{x,t- 


-Tl) 4 


- ^bpiAx,t) 


piA^^t) 


= apr{x,t) +apf 


{x,t) 


+ (1- 


b)piAx,t- n). 







We are interested in the diffusive behavior of the total probability p{x,t) = p^ + +P/.c that should 

be visible on large length and time scales. We therefore Fourier transform the master equations (|13|l and then Taylor 
expand all the coefficients up to first order in the frequency uj and up to second order in the wave number k. Thus 
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FIG. 8: One-dimensional random-walk model in a periodic sequence of gaseous and liquid regions. Five photon states with 
probabilities pj*, pj~, pp, p'f and pi_c are introduced. The subscript stands for the phase, in which the photon resides, and 
the superscript arrow for its direction of motion. Photons with angles Pi > /3*, i.e., photons undergoing total reflection in the 
liquid phase, are in the photon channelling state with probability pi^^. The possible transitions between the states are given by 
arrows. The transition probabilities and times are explained in the text. Note that the transitions from and into the photon 
channelling state are instantaneous. 



we obtain a set of linear equations: 
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(14) 



where /(r) = 1 — Ilot and = 1 — iiOTgi ± ifcAx ± 
ukTgiAx — k^Ax'^/2. Now, each probabihty in Eqs. itTH) 
obeys the same relation, e.g., detM p'^{u!,k) — where 
detM denotes the determinant of the coefficient matrix 
in (|14() . To leading order in oj and k, it reads detM c>c 
iu + Dk^. Thus, our model indeed reproduces diffusive 
behavior and gives a formula for the diffusion constant: 

^_ 2Ax^bT{l~ a){2- a) 

^ [Tg6(l -a)+Ti{a + b- 2ab)] [3aT 2a - AT + i] 

(15) 

(note that Tgi = (rg + t;)/2 was already used). 

We now make contact with our two-dimensional model 
foams. Photons in the four states p^, p^, pf* and 
correspond or lead to photons in the liquid phase whose 
angle with respect to the interface normal obeys /3 < 
P*, i.e., they do not exhibit total reflection. Thus we 
choose T — Tg^i where Tg^i is the average transmission 
probability introduced in section UlI Al 

In stationnary equilibrium, the flux of photons into 
and out of the photon channnelling state has to be equal: 
app -f app ~ bpi c- This implies a relation between a and 
b: 

a pi,c /j/'p(A)rf/3/ m 

- = — — ^ — = ^-TH. = 1 ■ (16) 

b Pi +Pi pmdPi % 



Here we have used the formalism of section lTlI Al to relate 
pi^c and pp to the probability p{Pi) and then have 

assumed an equal angular distribution of the photons, 
thus p{(3i) oc cos(/3i), as in section ITlI Al The quantity b 
denotes the probability for a photon to leave the photon 
channelling state. In our two-dimensional model foams 
this can only occur at a junction where three channels 
meet since in the channel itself reflections do not change 
(3 > f3* . The number of reflections necessary to reach 
a junction scales as d~^, we therefore assume b oc d/lo. 
Since the constant of proportionality is of the order of 
one, we take b — d/lo to compare the one-dimensional 
model to the simulation data. In Fig. |S1 e.g., variations 
in b will shift the region where the diffusion constant 
decreases. 

The lattice constant Ax determines the diffusion con- 
stant in the limit where the liquid region tends to zero 
(rf — > 0) . From Eq. lfT3|) one finds 

D{d^O)^^j^Axc. (17) 

In our model foams, Aa; corresponds to half of the typical 
bubble size. To compare the diffusion constant of Eq. 
(|15|l to our simulation results, we choose Ax and the 
times Ti and Tg as shown in Fig. |51 i.e.. Ax = \/?>Iq/2, 
Tg = ng{\/3lo — d)/c and t; = nid/c. As illustrated in 
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FIG. 9: The parameters Aa;, Tg and n of the one- dimensional 
photon-channelhng model in relation to a two-dimensional 
model foam: Ax = \/3lo/2, Tg = ng{\/3lo — d)/c and 
n = md/c. 

Fig. El and □ this gives an excellent agreement with our 
simulation data. So the one-dimensional model explains 
the main features of our approach to photon channelling 
very well. Absolutely crucial for the success of the model 
is the introduction of the photon-channelling state. 

V. CONCLUSIONS 

Photon channelling found in experiments by Gittings 
et al. is an interesting concept which increases our un- 
derstanding of the diffusive transport of light in foams. 
To contribute to its further study, we have simulated 
the transport of light in two-dimensional Voronoi foams 
based on ray optics. In perfect honeycomb structures, 
we find superdiffusive behavior which is connected to the 
motion of photons in the liquid channels. However, as 
soon as we introduce disorder into the foams, light trans- 
port becomes diffusive. This agrees with our finding that 
the diffusion constant D only exhibits a weak divergence 
for decreasing disorder parameter Sr. Furthermore, ac- 
cording to our simulations and in agreement with intu- 
ition, D decreases with increasing refractive index n/ and 
also with increasing channel-width-to-length ratio cI/Iq. 
We are able to model this dependence with the help of 
a one-dimensional random walk model but only when 
we introduce a photon-channelling state. Therefore, the 
importance of photon-channelling becomes evident. The 
decrease of D with increasing d/lo or liquid fraction e is in 
qualitative agreement with measurements by Vera et al. 

However, there is an important difference; whereas 
Vera et al. observe a divergence of the diffusion constant 
for e ^ 0, it assumes a finite value in our simulations 



and analytic model. The reason is that we consider inde- 
pendent transmission and reflection events at each single 
liquid-gas interface. If the film thickness becomes suffi- 
ciently thin, the effective reflectance and transmittance 
of thin films seems to be more appropriate, where inter- 
ference effects are taken into account ;.13| . They lead to a 
strong decrease of the reflectance for e — > and therefore 
to a divergence of D. 

Building on the work of Gittings et al. 0|, we are 
able to justify their observed relation between the pho- 
ton's path-length fraction / in the liquid phase and the 
liquid fraction e without any free parameters. Their fur- 
ther observation that photon channelling only occurs in 
the region 0.04 < e < 0.2 will help to improve our under- 
standing of light propagation in foams. For e < 0.04, the 
extension of the Plateau borders and films are probably 
too small that pure ray optics is applicable and therefore 
photon channelling breaks down. In this parameter re- 
gion, we suggest that the reflectance and transmittance 
of thin films which include interference effects, as already 
mentioned in the last paragraph, are the appropriate 
quantities to model the photons' random walk. How- 
ever, we have no understanding why photon channelling 
ceases for e > 0.2. Especially, our analytic considerations 
in Sec. IIII XI emplov very general arguments without re- 
lying on a concrete foam structure, so it is not evident 
why they should break down for larger volume fractions. 

Nevertheless, we show in this article that photon chan- 
nelling in liquid foams can be understood on the basis of 
ray optics through analytic considerations and simula- 
tions using simple model foams. It seems that details 
of the foam structure are unimportant. We are confident 
that the results presented in this article add another piece 
to the puzzle of explaining diffusive light transport in cel- 
lular structures such as aqueous foams. 
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